A dynamical chiral bag model 
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We study a dynamical chiral bag model, in which massless fermions are confined within an 
impenetrable but movable bag coupled to meson fields. The self-consistent motion of the bag 
is obtained by solving the equations of motion exactly assuming spherical symmetry. When the 
bag interacts with an external meson wave we find three different kinds of resonances: fermionic, 
geometric, and cr-resonances. We discuss the phenomenological implications of our results. 
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I. INTRODUCTION 

The MIT Bag model [|| and its chirally invariant ver- 
sions, such as the Chiral Bag model [|]|| and the Cloudy 
Bag model Q, continue to be useful tools in the study 
of the physics of the nucleon and other baryons. They 
have also been used extensively in the discussion of var- 
ious phenomena ranging from strange stars || to ultra- 
relativistic heavy-ion collisions || , even though these of- 
ten involve situations of high density/temperature where 
the applicability of the models is doubtful. 

In most of the bag model studies so far, because of 
its simplicity, a static spherical bag is assumed. The 
few notable exceptions, which allowed for the possibil- 
ity of a dynamical bag boundary, focused mainly on re- 
producing the correct phenomenological parity order of 
the low-lying states of the nucleon, although several ap- 
proximations and modifications to the theory had to be 
employed. For example, Rebbi and DeGrand M studied 
a bosonic bag and quantized the full system with per- 
turbation theory in the limit of small spherical oscilla- 
tions. The authors in Ref. |g|-|ic|| considered a fermionic 
bag with a surface tension, as well as the volume energy, 
and quantized only the motion of the bag boundary in 
the adiabatic approximation. Nogami and Tomio ]Tl| ] 
also quantized the motion of the boundary, but used the 
adiabatic approximation only for the mesons. Although 
these works gave a reasonable ordering of the low-lying 
states of the nucleon, the more fundamental question of 
whether it is consistent and feasible to use a dynamical 
bag to model hadrons was not addressed. That is the 
motivation of the present work. 

In a previous paper |T^| we proved that the original 
MIT bag model with massless quarks admits only one 
classical solution other than the static one, namely, a 
bag constantly expanding at the speed of light. We thus 
concluded that an additional field, such as the mesons in 
the Chiral Bag model, is needed to have a consistent and 
nontrivial dynamical bag model of hadrons. In this paper 
we implement a method that allows us to find the clas- 
sical solutions of a spherically symmetric chiral bag for 
any motion of the bag radius. In particular we look for 
the self-consistent solution of the full theory without any 



approximation, in which the motion of the bag surface is 
determined by the conservation of the total energy. We 
can thus study the full nonlinear features of the model. 

This paper is organized as follows. We first show the 
method we use to solve the problem. We then discuss the 
resonances found with a driven bag motion. In the third 
section the problem of the self-consistent surface motion 
is addressed, and we discuss the results obtained with dif- 
ferent incoming meson waves. We finally summarize our 
results and discuss their phenomenological implications. 
The appendix provides more details about the method of 
solution. 



II. METHOD OF SOLUTION 

The Lagrangian of the system we study is Q : 

£ = \ { [i (VV 1 ^ - (d^h^j) - B] 9 v (x)- 
1 vbia + ir •7r 75 )V'A s + d Al o-cPcr + <9 A1 7r-cP7rl , (1) 



u 

where 6y(x) is 1 inside the bag and outside and 



(2) 



A s being the surface delta-function. From it we derive 
the following Euler-Lagrange equations of motion: 



j^d^tp = inside the bag, 



(3) 



i^n^vb = — (a + if ■ Tr-y 5 )ip 

Jn 



on the bag surface, 
(4) 



(5) 



(6) 
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It is possible to look for spherically symmetric solutions 
of the above equations 0] by writing: 



-ia-rf(t,r) 



a = <j(t, r) 



7? = 7r(t, r)f 



(7) 



(8) 



(9) 



where u includes the spin and isospin parts and can be 
written as 



T,d> -| i,u>) 



(10) 



The arrows indicate the spins while u and d the up and 
down flavors of the quarks, and v satisfies 



(a + t)v = Q . 



(11) 



In Eqs. | and |Tl|, a are the three Pauli matrices and 
should not be confused with the field <j{t,r). Eq. [ll] en- 
sures that the RHS of Eqs. [| and |^ are spherically sym- 
metric and causes 7? to be radially directed, as in Eq. |^. 
This kind of solutions is hence called hedgehog solutions 

For a static bag an analytic solution is known j2|,|13|, 
which represents a stationary fermion field coupled at the 
surface of the bag to time-independent a and 7r-fields. 
Our goal is to find the hedgehog solution for any spheri- 
cally symmetric motion of the bag's surface. 

Substituting Eq. |^ for ip in Eq. || we obtain 



dt dr 



.dg df 2 
dt or r 



(12) 



(13) 



It is not difficult to verify |l2J that the general solution 
of Eqs. [12, 13 has the form: 
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g(t,r) = -[Q'(t-r)-Q'(t + r)} 



(14) 



d 2 ir d 2 n 29tt . 2 1 

rTT = — [g f + gf } S{R - r) . 



dt 2 dr 2 



r dr r* 



2U 



(17) 



For r 7^ R, we notice that a(t, r) obeys the equation of 
a free s-wave while ir(t, r), being the radial part of the 
vector 7?, satisfies the equation of a free p-wave. Hence 
we can look for a solution of the form 

a(t, r) = a in (t, r)6(R - r) + a ont (t, r) [1 - 6{R - r)] , 

(18) 



7r(t, r) = 7r in (t, r)9(R - r) + 7T out (t, r) [1 - 9(R - r)] , 

(19) 

where the fields inside and outside of the bag can be 
written accordingly as 

a in (t, r) = - [S in (< - r) - S fa (i + r)] + <r Q , ln (r) , (20) 



0- ou t(*,r) = - [E ou t-(t -r) - S ut+(i + r)] + cr . ou t(r) , 
r 

(21) 



7r in (t, r) = - <{ n( n (t - r) + n; n (t + r) + 



- [n in (t - r) - n in (t + r)] } + 7r ,i„ (r) , (22) 



Tout (*,r) = J |lCt-(* - r) + lCt + (* + r) + 
J [n out - (t - r) - n out+ (f + r)] J + 7r , out (r) , (23) 

where E out+ , S out _, £ in , n in , n out+ and n out _ are ar- 
bitrary functions. Notice that S ou t+ and S out _ are in 
general different functions as are also n out + and n out _. 
Here, the time-independent terms, <7o,in(?"); co.outM, 
7I'oin(' , ) and 7To out (f) ar e the static-bag solutions given 

by [HI 



f(t,r) 



Q'{t - r) + Q'(t + r) + 



[Q(t-r)-Q(t + r)] 



(15) 



where Q(z) is an arbitrary function. In spherical coordi- 
nates Eqs. H and ^ become 

d 2 G d 2 a 2 da 1 ,».,,,„ , f . 

W-w-;-b7 = -2U [99 - ff]s{R - r) ' (16) 



O"0,m(f) = go , 

ffo.out (r) = go + olRq 



1 

#0 



7T0,in(f) 





-— r 
3 



7T0,outW = ■ 

3 r* 

Substituting Eqs. |TJ and [TJ in Eqs. |TJ and |TJ and re- 
quiring the continuity of a(t, r) and ir(t, r) at r = i?, we 
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finally obtain two relations that can be viewed as bound- 
ary conditions for the fields g, f, 

&'m j ^out 7 " in 

and 7T ut: 



/,» [ da ' in _ da ° ut \ _|_ ( d<J ' m _ dcout ^ _ 



dt 



dt ) 

~{g*g-rf) 



<9r 



dir„ 



V at 



(9f 



dr 



dr J 
at r = R , 



dir a 



(24) 



^ (<?*/ + <?/*) 



9r 

at r = R 



(25) 



From Eq. || we can express a(t, R) and 7r(t, R) in 
terms of g(t,R) and f(t,R) (see Appendix) and use 
Eqs. |24| and ^5] to find g(t,r), f(t,r), i.e. Q Ie (z), and 
Qim(z) (see Eqs. [l4|, |l5|) with the null-lines method 

EMI- 

From the point of view of the null-lines method the 
unknowns in Eqs. [24| and as long as \R\ < 1 |fj"5| , are 
Q rc (t+R), Q im (t+R), £ in (*+#), Il ia (t+R), Eout- (*-#), 
and n ou t-(t — -R). Furthermore, the latter four are fixed 
once Q TC (t + R) and Q lm (t + R) are known, using 

S in (t + R) = E in (t-R) + Rg -Ra(t,R) , (26) 



E out _(* - i?) = E out+ (t + R) + aRg 



R(g + aR )+Ra(t,R) , 



(27) 



n? n (t + i?) = 



1? 



[n ln (t + 
i?) 



fl) - n in (t - i?)] + 

, /3 



-If 



Rir(t,R) , (28) 



n; ut _(f-i?) 



1 

R 



[n out+ {t + R)-n out -(t-R)} - 

(3Rl 



n' 



out+ (t + R) + ^+R7r(t,R). (29) 



Here for convenience we have not written explicitly the 
dependence on Q IC (t + R) and Q lm (t + R) which are hid- 
den in a(t, R) and ir(t, R). 

We still need to solve Eqs. ^8] and ^9|. This can be 
done numerically either by simply replacing Tl[ n {t + R) 
and TL' out _(t — R) with their finite-difference counterparts 
or by integrating between z — dz and z, where z = t + R 
and z = t — R respectively for the first and second equa- 
tions, and by approximating all quantities other than 
IT in (t + R) and II out _(f — R) as being constant in this 
infinitesimal interval. The numerical results turn out to 
be slightly more accurate with the second method. To 
solve Eqs. and |2^ we used a fourth-order Runge-Kutta 
algorithm. 



III. RESONANCES WITH A DRIVEN BAG 
MOTION 



With the method discussed in the previous section we 
first computed the solution for a static bag and then for 
a slowly moving one. We verified that the norm of the 
fcrmion field is conserved and that our numerical method 
is accurate up to the second derivative of Q(z) for a bag 
of initial radius i?o = 1 f m and /„■ = 1 fm^ 1 . With these 
parameters the static chiral bag is similar to the MIT 
bag, with an almost constant o~(r) and a very small ir(r). 
In Fig. |l| we plotted the second derivative of Q(z) in or- 
der to show the quality of the numerical solution. All the 
results presented below are obtained with R = 1 fm and 
= 1 fm -1 , which is a representative set of parameters 
for showing the qualitative features of a dynamical chiral 
bag model. 

Since we are particularly interested in the behaviour 
of the fields under the effect of the motion of the bound- 
ary, we first study the chiral bag with an imposed sur- 
face motion. Subjecting the bag boundary to a sinu- 
soidal motion, R(t) = Rq + e[cos(ut) — 1], we found 
three different kinds of resonances: i) the fermionic res- 
onances, which are excited when the oscillation frequen- 
cies are close to the difference between two static-bag 
eigen-energies, v ~ E n — Ek', ii) the geometric a reso- 
nances, for v ~ mr/Ro; iii) the parametric a resonances, 
for v ~ (2n + l)ir/(2Ro), where n is an integer. 

The origin of the fermionic resonances at v ~ E n — Ek 
is similar to those found for a Schrodinger particle in an 
oscillating cavity The difference here is that the 

fermions cannot really be excited to the upper static- 
cavity level because the upper level is associated with 
different static pion fields which cannot be produced by 
the boundary motion. However, since with our choice 
of parameters <r(r) and 7r(r) change little for different 
static solutions, the system still gets excited for oscil- 
lation frequencies close to the static energy gaps. The 
smaller f v is, these resonance frequencies deviate more 
from E n — Ek- As an example, we show for v = 3.4/i?o — 
E-2 — Ei the time-dependence of the energies of the 
fermion and the meson fields in Fig. ^a, b and c respec- 
tively. It is interesting to note that neither the a nor the 
7r-fields gain considerable energy. 

For v = nir / Rq we found resonances involving the er- 
field. As can be seen in Fig. || the energy of the tr-field 
increases remarkably, while the energies of the fermion 
and the 7r-field change little. These resonances may be 
considered geometric, since the resonance frequencies are 
related to the time it takes for the wave components of 
a(t, r), i.e., Ti- m (t — r) and Tt- m (t + r), to travel from the 
boundary of the bag to its center and back again. It 
has been shown [p]] that p-waves in an oscillating spher- 
ical cavity also manifest resonances at v — mr/Ro, an d 
so it is somewhat surprising that here the energies of the 
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fermion and 7r-fields are little affected. The strongly non- 
linear interaction at the bag boundary seems to damp out 
the resonant evolution. The energy of the fermions actu- 
ally shows some resonant behaviour, but this is probably 
mainly due to the fact that the driving frequency is close 
to E 2 - E x . 

The third kind of resonances we found is a peculiar 
feature of the system under analysis. As we can see 
from Fig. |[ it involves mainly the er-field. Note that 
Sj n (z) (Fig- has an almost periodic dependence and 
the period is about half that of the oscillating bag. In 
other words the bag surface, oscillating at frequencies 
v = (2n + l)ir/ (2Rq), acts as a source for a cr-field with 
frequencies (2n + 1)tt/Rq. Such frequencies are obviously 
resonant with the cavity so that the a field is resonantly- 
enhanced. The occurence of this kind of resonances can 
be understood in the following way. Since the fermion 
field is not excited we can approximate its value at the 
bag boundary with its static-bag expression 



g(t,R(t)) = Nex V {-iEt)j {ER(t)) , 



(30) 



f(t, R(t)) = -iVcxp {-iEt)h(ER(t)) . (31) 
3|, we obtain after some manipu- 



From Eqs. g ||, | 
lations 



fl(t,R) + f? m (t,R)-gl(t,R) 



Sim 



(t,R) 



fl{t, R) + fUt, R) + gl(t, R) + g* m {t, R) 



a(t,R) = -2f„ 



(32) 



frejt, R)g IO (t, R) + f im (t, R)g im (t, R) 

fl(t, R) + /&,(*, R) + <&(*, R) + 9L(t, R) ' 



(33) 



One can see that the dependence on exp (—iEt) cancels 
out in our approximation and the whole expressions be- 
come periodic functions with period T = 2ti/v. The 
Fourier expansion of such functions involves all multiple 
frequencies of v and, in the case of v — (2n + l)ir/(2Ro), 
its even multiples are also integral multiples of it/Rq. It 
is then evident that the expressions for cr(i?) and 7r(i?) 
contain terms in resonance with the cavity. However, it is 
surprising that for oscillation amplitudes e > 0.005i?o the 
frequency 2v for a[R) becomes the dominating one even 
before the first bag oscillation is completed. In Fig. |] we 
can see clearly how the amplitude of Si n (z) increases with 
each bag oscillation. Although the expression for tt(R) 
also contains terms with frequency being integral multi- 
ples of it/Rq-, we observe no resonant behaviour for the 
7r-field, which is consistent with the previous observation 
that this field is not excited for v = nir/Ro- 



IV. SELF-CONSISTENT SURFACE MOTION 



A. Equation for the radius and non-resonant 
interaction 



Our main interest in this work is to study the be- 
haviour of the fields and the bag's surface when perturbed 
from their static-bag states by, for example, an incoming 
pion wave. To this end we need to find the self-consistent 
dynamics of the bag surface and fields. 

We notice that the velocity of points on the bag sur- 
face does not appear in the Lagrangian, Eq. [|, and hence 
we cannot derive from it an equation of motion for the 
bag radius Jj|H • The motion of the bag however is con- 
strained by the conservation of the total energy. Let us 
consider the energy-momentum tensor 

T" v = -g^C + l - [^d u ^j - d v ^^\ 9 V + 

d^od v o + d"n ■ d u n . (34) 

The conservation of energy and momentum requires 
d^T^ — 0, from which, by using the equations of motion 
Eqs. g, |, |, | and after some algebra, we derive 

Bn v = ——d u [tp((T + if ■ 7?75)-0] on the surface . 

(35) 

For spherically symmetric solutions n v = (i?, f), and the 
above equation can be written as 



BR = 1 ^ [${<r + ir ■ iVya)^] j 



r=R 



B=—l— \jj(<T 

2f n 1 dr L 



(36) 



(37) 



r=R 



In the static case Eq. |3^ is identically satisfied because 
R = and ip{a + if ■ f^)ip is time independent, and 
we could use Eq. [37] to derive B. However, the RHS of 
Eq. H?] is an ambiguous expression, because it involves 
the derivatives of a and 7r at the boundary which are dis- 
continuous. To overcome this difficulty we use the fact 
that T^ v can also be written as [0 

T^ = T^6v + TZ(l-e v ) , (38) 

because the surface term is zero along the trajectories of 
motion. Since d^T^ = and d^T^t = 5 the conserva- 
tion condition for energy and momentum becomes 

on the surface . (39) 

Again using the equations of motion we obtain 



"/* in — 'V I out 
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with 



j^d" (er av + if • 7? av 7 5 ) 1p] 

-<V<9'V in <9 , V ou t + n l _ l d tl a ou td v 'a in 
-n^ir in d v K out + n fl d tl Tr out d' / TT in = , 

- 1 r ^ 

Cav = g [0"ia + C ou t j ! 



TTav = 2 fan + ''"out) 



(40) 



(41) 



(42) 



(9 p (Ji n ) 2 + ((9 p 7r in ) 2 - (<9pCr out ) 2 - (9 p 7f out ) 2 



(43) 



and all the functions in the above expressions are evalu- 
ated at the surface of the bag. 

Eq. ^ is well defined, and since the spatial part of n v 
in the case of spherical symmetry is f , it can be used to 
calculate the bag constant B as 

19- 

B = D(t) + TTTjr ( CT av + IT ■ 7?av75) lf\ 

nu da out da m 

-n,xO^CT in — 1- Tl^U O out 



dr 



dr 

au dir out au diri 



dr 



'dr 



(44) 



where the static-bag solution has to be used. With the 
static hedgehog solution the terms involving products of 
fields inside and outside the bag actually cancel out each 
other. Such a value of B ensures that the continuity 
equation for the linear momentum is satisfied. Choosing 
v = from Eq. ^o| we can derive an equation for R: 



-1 



' 1 

]r= B—W){2f,dl^^ 
da hl da out 
dr dt 
dir in dir out 
dr dt 



■ IT ■ 7? av 7 5 ) 1p] 

gqWt d<7 in _ 

dr dt 
dir out dir in 

dr dt 



(45) 



We use the above expression to compute the motion of 
the bag's surface, which conserves the total energy and 
momentum. 

It is important to notice that for spherically symmetric 
solutions the total linear momentum is conserved regard- 
less of the value of B, because the associated current is 
radial and the vector sum always gives a zero total mo- 
mentum. This guarantees the conservation of the total 
momentum also for a non-static bag surface, because in 



such a case the RHS of Eq. |4| is not constant and hence 
the equation is not satisfied. 

The first question we want to address is whether the 
static hedgehog solution is stable with respect to a small 
perturbation or it is just a special field configuration per- 
mitted only with a static boundary. If the static hedge- 
hog models a hadron state one would like it to be little 
affected by a small non-resonant perturbation. We there- 
fore considered an incoming wavepacket incident on the 
bag, and we computed the evolution of the system. Wc 
used both a 7r-field and a a- field as the incoming packet, 
with the following form 



C ou t(*0) r ) 

TTout {to, r) 



= A[e-^ r - R ^ -I] 3 - 

sin (v(t + r)) e -«('--Ro) r > i? , 

(46) 



which, at r = Rq, is zero up to the third derivative, in 
order to avoid discontinuities at the instant of the colli- 
sion. The motion of the bag surface depends on the bag 
constant B. We use B = 1 fm to demonstrate the 
qualitative features of the system. 

In Fig. H we show the energy of the fields inside and 
outside the bag versus time for small A and v. In all 
the computations we used (3 — 1/Rq. Both for a 7r-wave 
and a cr-wave the bag is hardly changed, and after the 
interaction the velocity of the surface goes back to zero 
gradually as expected. Part of the pion field is reflected 
back at the surface of the bag, while the other part after 
penetrating the bag goes back out towards infinity. It is 
interesting to observe that a cr-wave has almost no effect 
at all on the bag and nearly does not enter it. We verified 
that the static hedgehog solution remains little affected 
also for larger A and higher frequencies v, thus validating 
its use as a stationary state of a hadron. 



B. Resonances 

We next consider whether the resonances found in 
the case of a driven bag motion still occur for the self- 
consistent motion caused by an incoming wave of ap- 
propriate frequency. This is a non-trivial question be- 
cause the nonlinear relation between the fields and the 
motion of the boundary, as expressed in Eq. ^|, might in 
principle destroy any phase coherence on which a reso- 
nance is built up. We hence performed our computation 
with incoming 7r-waves given by Eq. ^ with a = and 
V = nir/R , v = (E n - E k ), and v = (2n + l)vr/(2i? ). 
Again due to numerical limitations we had to use small 
values of A. 

In Fig. U we plot the energy of the fields inside the bag 
vs. time for v = tt/Rq and in Fig. [?] for v = E2 — E\. 
We can see that the resonant behaviour is still present. 
From the point of view of the energy gained by the bag 
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the two resonances merge and appear as a broad reso- 
nance peaked at v = — E\. With a closer look how- 
ever one can observe two different physical phenomena. 
For v close to 7r/i? the c-field is excited and the bag 
expands while the fermion field gets only a small con- 
tribution from the second static hedgehog state. For v 
close to E2 ~ Ei the er-field is slightly out of resonance 
while the fermion field tends to be excited towards the 
second static hedgehog state and the volume of the bag 
decreases. 

In the case of a wave with v — tv/(2Rq) we also have 
a clear resonance that involves the a- field (Fig. ||). At 
this frequency the fermion field is completely out of reso- 
nance. Overall the bag's energy increases not only as the 
pion energy but also in the form of volume energy due 
to the expansion of the bag. The excitation mechanism 
for the cr-field is the same as explained in the previous 
section by means of Eqs. [32| and [33| 

It is very interesting to notice that the expansion of 
the bag is related to the excitation of the a-field and 
not directly to the fermions. The increase of energy due 
to a larger bag radius is the classical counterpart of the 
breathing modes proposed by other authors j7|-[iT| to ex- 
plain certain radial excitations of the baryons such as the 
Roper resonance. These authors propose that such reso- 
nances are excitations of the collective degrees of freedom 
of the bag, represented in the models they considered 
by the surface coordinates, while the quarks essentially 
remain in the ground state. In the Chiral Bag model, 
besides the bag's radius, the pions too describe collec- 
tive degrees of freedom, and hence our results strongly 
support the above scenario. 

For larger odd multiples of tt/(2Rq) the resonant be- 
haviour is much attenuated. We believe that this is due 
to the fact that at higher frequencies and with a self- 
consistent bag motion an approximation as the one shown 
111 Eqs. |^ and H is no longer acceptable. The perturba- 
tion still appears to be resonant for the a field, but its 
energy increases very slowly. 

Another remarkable property of the chiral bag is that 
it shows a realistic behaviour in a scattering process. We 
have already mentioned the dynamics for the scattering 
with a non-resonant pion wavepacket, but it is interesting 
to examine how the bag releases its energy after being ex- 
cited by a resonant wavepacket. In Figs. ^ it can be seen 
how, after the incoming wavepacket has been scattered 
away, the bag remains in its excited state for some time 
before starting to slowly release the energy gained, hence 
looking like a stable particle. Also in this case we can 
see that the 7r-wave inside the bag is not excited and its 
temporary increase in energy is simply due to the part 
of the incoming wave that enters the bag before being 
reflected away. 

Since the chiral bag is generally used to describe 
baryons, we have also performed the previous calcula- 
tions with three quarks inside the bag in order to ver- 



ify that our findings hold in this case too. Not having 
quantized the theory, we have to consider already in the 
Lagrangian three distinct fermion fields. Moreover, the 
coupling with the pions causes the solution to differ from 
the one-quark bag. In fact, while each quark has to sat- 
isfy the same equations Eqs. || and [|, the equations for 
the pions change because the RHS of Eqs. || and |^ must 
be multiplied by a factor 3. Such a difference produces 
different eigenvalues for the static cavity solution and 
yields the following equation for the motion of the bag's 
surface 




We have verified that the features found with only one 
quark remain with three quarks. 

For smaller values of B or bigger amplitudes A we have 
been able to obtain fairly accurate solutions for short 
times (t < 10 fm/c) and, for resonant perturbations, we 
observe a much stronger and faster excitation process, 
which could probably model realistic energy levels. 

V. SIMILARITIES WITH THE SOLITON 

We point out here the very interesting similarities be- 
tween the hedgehog solution and the soliton solution of 
a non-linear field theory. 

We can define a classical soliton as any spatially con- 
fined and nondispersive solution of a classical field theory 
pit . In order to have soliton solutions it is necessary to 
have some nonlinear couplings among the fields. The 
MIT bag model with only fermions inside the bag does 
not have nonlinear couplings. In fact, as we proved in 
a previous paper p"2| , it admits bag-like solutions, but 
these are unstable with respect to perturbations of the 
bag surface. In the chiral bag model the quark-pion cou- 
pling, although it is linear, introduces a nonlinear self- 
coupling for the fermion field through the boundary con- 
ditions, as apparent from Eqs. [24|, |25|. We have seen that 
the hedgehog solution is indeed stable with respect to 
perturbations of the bag surface. 

For a boson field the various nonlinear couplings can be 
characterized by a dimcnsionless coupling constant g. If 
g = 0, the theory is linear and there is no soliton solution. 
However if <?, however small, is different from zero, the 
theory admits soliton solutions. In the limit g — > the 
soliton solution grows to infinity. This is remarkably sim- 
ilar to what happens in the dynamical chiral bag model. 
The dimensionless coupling constant is in this case 7//^, 
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with 7 an arbitrary constant with dimension of L _1 . If 
we set 7//tt = already in the Lagrangian, we have the 
MIT bag model, and no stable, spatially confined solu- 
tion exists. If we take the limit 7//^ (/„• — * 00), we 
still have stable hedgehog solutions, but the field a goes 
to infinity, so that a/ / w — > 1. 

Such similarities seem more than a coincidence, espe- 
cially if we consider the bag models as simplifications 
of more general models. In fact it has been shown JToj ] 
that a chiral model, similar to the Skyrme Lagrangian, 
can automatically produce bag-like solutions. From this 
point of view one may put some features of the bag- 
like solutions already in the Lagrangian, thus obtaining 
a bag model. The fact that the MIT bag model does 
not admit a stable bag-like solution may be viewed as 
due to an oversimplification, having completely neglected 
the quark-pion interaction, while in the chiral bag model 
such interaction is maintained at least at the surface of 
the bag. 

VI. CONCLUSION 

We have shown that the chiral bag described by the 
Lagrangian Eq. [I] is stable in the sense that a bag-like 
solution exists even if the static bag is perturbed either 
by arbitrary radial motions of its boundary or by its in- 
teraction with a meson wave. This is in contrast with 
the purely fermionic MIT bag which has been proved to 
be unstable 

Computing the solution for a bag perturbed by a non- 
resonant meson wave, we found that it remains close to 
the static hedgehog solution and that, after the incoming 
wave is scattered away, it returns to the static hedgehog. 
Such a result validates the use of the static hedgehog as 
a stationary state of a hadron. 

We also examined the existence of resonant perturba- 
tions, and we discovered three kinds of resonances which 
occur when the bag interacts with an incoming 7r-wave. 
When the frequency of the incoming wave is close to an 
energy gap, v ~ E n — Ek, the fermion field is in reso- 
nance. The cr-field is also excited since v is close to an 
integral multiple of tt/Rq. The fermion field tends to go 
to an upper static hedgehog level, but these resonances 
are not simply transitions from a lower hedgehog state to 
an upper one, because in that case there would be only 
a static meson field in the final state, while here we have 
remarkable energy contribution from a non-static er-field. 
At v = nir/Ro the fermion field is little excited while the 
(T-ficld is in resonance and the bag expands. The third 
type of resonance occurs when v equals odd multiples of 
7t/(2Rq), which is a consequence of the linear boundary 
condition Eq. ^. In this case the fermion field is not ex- 
cited and the increase in energy comes from the cr-field 
and the increase in the volume energy associated with 
the expansion of the bag. The occurence of a resonance 



at v = 7r/(2i?o)j which is much smaller than the energy 
gap between the first and second static hedgehog states, 
gives support to the description of the Roper resonance 
as a radial excitation of the collective degrees of freedom. 
In particular the expansion of the bag without fermion 
excitation is the classical analog of the breathing modes 
proposed by other authors J7| pr[ . Our results, however, 
show that the expansion of the bag is strictly related to 
the excitation of the cr-field and in this sense warn against 
the use of the adiabatic approximation. 

Compared to previous studies of dynamical bag mod- 
els, our approach has two main advantages: we solve the 
equations of motion without approximations and we show 
the dynamics of the resonances. 

In the present work we have considered hedgehog con- 
figurations in which the quarks are neither in a flavor 
eigenstate nor in a J-eigenstate; they do not represent 
any known hadrons. However the occurence of the reso- 
nances we found does not depend on the flavor of the 
quarks, and since they are caused by spherical waves 
the angular momenta are not changed. Moreover they 
are not related to specific features of the hedgehog so- 
lution. In fact the v — E n — Ek resonance is a gen- 
eral feature of discrete-level systems. The resonances at 
v = nn/Ro are geometrical resonances related to the fact 
that the mesons in this model are massless, and the ones 
at v = (2n + l)n/(2R ) are related to the linear bound- 
ary condition Eq. 0. Therefore, we conjecture that such 
resonances should occur also for more realistic solutions, 
i.e. with definite flavor and J. 

We thank the support of a Hong Kong Research Grants 
Council grant CUHK 4189/97P and a Chinese University 
Direct Grant (Grant No. 2060193). 

VII. APPENDIX 

For a spherical bag n M can be written as (i?, — f) and, 
as shown for example in jl2j , Eq. |] in radial coordinates 
becomes 

iitg(t, R) - f(t, R) = ±- [<r(t, R)g(t, R) - n(t, R)f(t, R)\ , 

J TT 

(48) 

- iitf(t, R) - g(t, R) = ^~ [tr(t, R)f(t, R) + n(t, R)g(t, R)] . 

(49) 

Equating separately the real and imaginary parts of the 
two equations we obtain the following four equations 

-Rg im (t,R)- f IO (t,R) = 

(50) 

f [a{t, R)g ie (t, R) - n(t, R)f m (t, R)} , 
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Rfun(t,R)- g lc (t,R) = 
± [a{t, R)f ie (t, R) + w(t, R)g re {t, R)] 
Rg Ie (t,R)- f im (t,R) = 
j- [a(t, R)g im (t, R) - 7r(t, R)f im {t, R)} 
-Rf Ie (t,R)-g im {t,R) = 
f [o-(i, fl)/ im (t, i?) + 7r(t, R)g im (t, R)} 



(51) 



(52) 



(53) 



Equating Eq. p4 with Eq. foB and Eq. ml with Eq. 57 we 



/»(*, fl)/im(t, R)]R I (gl(t, R) + &{t, R)) , 



(55) 



w(t,R) = U f? m (t,R)~gf m (t,R) - [f rc (t,R)g im (t,R)+ 



(56) 



7T(t, i?) = /» \ fUt, R) ~ <&(t, R) + [fre(t, R) 9im (t, R) + 



obtain two non-linear equations for the real and imagi- 
nary parts of / and g, 



2f im (t,R)g im (t,R) - [g IC (t,R)g im (t,R) 



/«,(*, R)f im (t, R)] R> ■ [ fl *,(t, R) + f? c (t, R)] 



(58) 



At this point we have four equations and four un- 
knowns, i.e. Q re (t + R), Qim{t + R), cr(t, R), and n(t, R). 
This fact seems to make our requirement, that <r(t, r) and 
n(t, r) be continuous at r = R, redundant hence making 
the whole problem inconsistent. In fact if we could de- 
rive from Eqs. |5(] , ^l] , [5^ and |5^ all the four functions 
mentioned above, then Eqs. [l6| and [Tt] would require cr ou t 
and 7r ou t to be singular at r = R, and even this would 
not guarantee the existence of a solution in general. How- 
ever, it turns out that Eqs. |5(] , |5l] , [32] and |5^ arc not 
independent and we need to impose some condition on 
the functions in order to have a unique solution. We have 
verified this in two ways, as explained below. 

Solving Eqs. §I|, || and || for a(t, R) and ir{t, R) 
we obtain two different expressions for each function 

a{t, R)=f^\- 2f im {t, R)g im {t, R) + [g ie (t, R)ton(t, R) 



(54) 



a(t, R) = f„\ - 2/ re (t, R)g ie (t, R) - [g re (t, R)g im (t, R) - 



2f re (t, R)g ie (t, R) + [g IC (t, R)g im (t, R) 



/»(*, R)fUt, R)]r\- [gL(t, R) + fL& R )] 



f? m (t, R) - gf m (t, R) - [f ie (t, R) 9im {t, R) + 



gre(t, R)f im (t, R)]R\- [gUt, R) + /£(*, R)] 



(59) 



fl(t, R) - gl(t, R) + [f Ie (t, R)g im (t, R) + 



g ie (t, R)f im (t, R)]R \- [g? m (t, R) + R)] 



9re(t, R)fun(t, R)]r\/ (gfjt, R) + tfjt, R)) , 



ffre(t, R)fim(t, R)]R\ / (fl£,(t, R) + fl{t, R)) 



The problem is evidently extremely difficult to han- 
dle analytically, and so we used a numerical ap- 
proach. Substituting Q' with a finite incremental ratio, 
Eqs. |58| and ^9] become two non- linear algebraic equations 
for Q Te (t + R) and Q lm {t + R). Solving numerically the 
algebraic equations, we have found that a whole region 
exists, in the Q rc — Q lm plane around Q(t + R — dz), 
in which the algebraic equations are satisfied, thus indi- 
cating that the solution is not unique. Since this is not 
a rigorous proof, we need to cross-check our finding by 
imposing the continuity of a(t,r) and tt(£, r) on the sur- 
face of the bag and by verifying whether the solutions, 
found without using Eqs. ^8] and satisfy all the four 
Eqs. |(J [51], H and |§ 

In order to do this we have solved numerically 
Eqs.||Jand|2^for Q le (t+R) and Qi m (t + R), as discussed 
in the main part of the paper, with a(t,R) and ir(t,R) 
given by anyone of the above expressions or a combina- 
tion of them. We thereby verified that Eqs. [SO], ||l], [32], |53| 
are automatically satisfied. 

In the case of R = we can analytically prove that 
Eqs. |5(], 51, 52 and 53 are not independent. In fact, 

P(t)g(r), fre = 



(57) looking for solutions of the form g r 
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p {t)f( r ), 9im = S(t)g(r), f im = S{t)f(r), we easily ver- 
ify that Eqs. ^9] are automatically satisfied leaving 
P(t) and S(t) undetermined, ft is reasonable to think 
that as R slowly departs from zero the four equations 
still remain not independent, though we are not able to 
provide a rigorous proof for it. 
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FIG. 1. Time evolution of the real and imaginary parts 
of Q"(z) for driven surface oscillations, with v — 1/Rq and 
e = 0.05i?o. 
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FIG. 2. Time evolution of the energy for driven surface os- 
cillations at the v = tt/Ro (dashed lines) and v = E2 — E\ 
(solid lines) resonances with e = 0.01-Ro- a) Energy of the 
fermion field, b) Energy of the 7r-field. c) Energy of the 
(7-field. 
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FIG. 3. Time evolution of the energy for driven surface os- 
cillations at the v = 7r/(2i? av ) resonance with e = 0.08i?o- 

(-Rav = Ro + e) 
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FIG. 4. The function Eia(z) for driven surface oscillations 
at the v = 7r/(2i? av ) resonance with e = 0.08i?o- The function 
clearly contains a periodic contribution with T ~ 2/Rq. 
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FIG. 5. Time evolution of the energy in the case of a colli- 
sion with a wavepacket with v = 1/Ro, a — 0.2/Ro, 7r-wave 
with A = 0.05/i?o (dashed line) or cr-wave with A — 0.1/ Ro 
(solid line) (see Eq. kq): a) volume energy, b) fermion energy, 
c) 7r-field energy inside the bag, d) a-field energy inside the 
bag, e) 7r-field energy outside the bag, and f) cr-field energy 
outside the bag. 
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FIG. 6. Time evolution of the energy for a resonance with 
an incoming v — tv/Rq, A — 0.005/i?o, a = 7r-wave (see 
Eq. ^6j), showing a) fermion energy plus volume energy, b) 
7r-field energy inside the bag, and c) cr-field energy inside the 
bag. 



M 



pi 

b 

w 




FIG. 7. Same as Fig. 6, but for v = Ei-Ex. 
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FIG. 8. Same as Fig. 6, but for v = tt/(2R ) and 
A = 0.02/i? - 



+ 



6.2 

0.15 

0.10 

0.05 - 

0.00 
0.0015 - 
0.0010 - J 



0.0005 
0.0000 





b) 



c) 



40 

t/R 



FIG. 9. Same as Fig. 6, but for = 6.55/i? ~ E 3 - E u 
A = 0.005/i? and a = 0.2/ R . 
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